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We derive an analytic expression for the 0(a) radiative correction to the u e (u e ) 
spectrum in allowed /3-decay. The m — > limit is convergent and leads to a very 
simple result (m is the electron mass). This is in sharp contrast to the correction to 
the e~(e + ) spectrum, that diverges as m — > 0, an important difference that we explain 
on theoretical grounds. After discussing some of their general properties, we use the 
corrections to the v e and e~ spectra to derive the corresponding correction to the 
conversion from the e~ spectrum to the v e spectrum, a relation that is of considerable 
interest for reactor studies of neutrino oscillations. 
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1 Introduction 



The 0(a) radiative correction (R. C.) to the e~(e + ) spectrum has played an important role in 
precision studies of /3-decay. It is described by a R. C. (a/27i)g(E, E ), where E and E are the 
energy and end-point energy of the e~(e + ) [T|2] . Here we have in mind the R. C. not contained 
in Fermi's Coulomb function F(Z,E), which is separately included in the theory of /3-decay. 
The effect of (a/2ir)g(E, Eq) on the spectrum shape is significant, particularly for large values 
of Eq [T1I3] . It also gives important contributions to the theoretical calculations of the Ft- values 
of superallowed (0 + — > + ) Fermi transitions jl] and the neutron lifetime [5]. We recall that the 
superallowed Fermi transitions play a crucial role in the precise verification of the unitarity of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. In fact, the current analysis, which includes large 
short-distance effects jl], the R. C. (a/2ir)g(E, E ), and a number of more recent refinements p]6] , 
leads to the conclusion [7]: 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.99995(61) , (1) 

an impressive test of the three-generation Standard Theory at the level of its quantum corrections. 
It is also interesting to note that the overall electroweak corrections to Eq. ([1]) are of 0(4%), i. e. 
66 times larger than the 0.061 % error. 

The correction (a/27r)g, in conjunction with the short- distance effects and refinements dis- 
cussed in Refs. [^|5f6], has also been recently applied to the calculation of the muon capture rate 
in the singlet state of H [8]. 

The aim of the present paper is to derive the corresponding R. C. to the v e {y^) spectrum in 
allowed /3-decay which, as explained below, is of considerable current interest. 

An important difference with the R. C. to the e~(e + ) spectrum involves the inner bremsstrahlung 
(I. B.) contributions. In fact, the integration over the e~(e + ) momenta introduces a complicated 
dependence on the radiated photon energy that hinders the subsequent integration over the photon 
momenta. In Section [2] we develop and implement a convenient strategy to deal with this problem. 
The combination of the I. B. and virtual corrections leads then to the R. C. (a / 27r) h(E, E ) to 
the v e [y^) spectrum, which we present in analytic form. Here E = Eq — K, where K is the P e (f e ) 
energy. We find that in the m — > limit (m is the electron mass), Ii(E,Eq) converges and has a 
very simple form. This is in sharp contrast with g(E, Eq), that diverges as m — > 0, an important 
difference that we explain on theoretical grounds. 

In Section [3] we apply the R. C. (a/2ir)g(E, Eq) and (a/2ir)h(E, Eq) to the conversion from 
the e~ spectrum to the v e spectrum, a relation that is of considerable current interest for reactor 
studies of neutrino oscillations [9"|10|lllj . This leads to the analytic expression for the 0(a) R. C. 
to that important relation, which we illustrate with numerical examples. 

2 Radiative Correction to the v e (v e ) Spectrum 

As in Ref. [T], we first consider the basic process of neutron /3-decay, n — > p + e~ + u e , and 
neglect terms of 0(aE/M), where M is the nucleon mass. 

The 0(a) virtual R. C, as well as the zeroth order transition probability, are obtained from 
those in the e~ spectrum by simply replacing E — > E,p — > p, (5 — > (3, where E,p, and m are the 
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e energy, momentum and mass, (3 = p/E, K is the v e energy and E,p, and /3 are defined by 

E = E -K; p=(^E 2 -m 2 y /2 ; /? = p/E . (2) 

The reason for this result is that, in deriving the v e spectrum, one integrates over the e~ phase 
space factor d 3 p/2E rather than the u^s d 3 K/2K. In these two contributions the conservation 
of energy is expressed by S(E — K — E) and the phase space integrations lead to the relation 
explained above. 

We now turn our attention to the inner bremsstrahlung (I. B.) contributions, which present 
more subtle and theoretically interesting differences with the e~ spectrum. 

Calling e 1 the radiated photon energy, the I. B. amplitude involves two terms: the first one 
behaves as l/e 7 as e 7 —> 0; the second one is of 0(e®). The square of the absolute value of the first 
term, combined with the d 3 k/2e 1 photon phase space factor (k is the photon momentum) leads 
to the infrared (I. R.) divergent contribution. The interference of the two terms and the square of 
the absolute value of the second one are I. R. convergent and are evaluated separately. 

As in [lj, we regularize the I. R. divergence by introducing an infinitesimal photon mass A min , 
so that e 1 = (k 2 + A^ in ) 1//2 . After integrating over the proton and electron momenta, as well as 
the v e direction, the I. R. contribution is proportional to J^ m (d 3 k/2e 3 )f(e 1} v ), where 

/(e 7l v) = pE { (1/pv) In [(1 + p v )/(l - 0v)l - (1 - /? 2 )/(l - /?V) - 1} , (3) 

v = k/ej, k m = E — m is the maximum photon energy for given K, and it is understood that 

E = E - e 7 , (4) 

which follows from energy conservation. 

The £ 7 variable in f(s 1 , v) denotes the e 1 dependence contained in E via Eq. (BJ and, therefore, 
also in p, (3, and Eq. (jHJ). The complicated dependence of f(e 1 ,v) on e 7 and v hinders the 
subsequent integration over the photon momenta. To deal with this problem, a convenient strategy 
is to split f(e 1 ,v) according to 

f{e^v) = f(0,v) + [f(e 1 ,v)-f{0,v)} . (5) 

The I. R. divergence arises from /(0, v), a function that only involves the photon momenta through 
its dependence on v. Changing the integration variable from k to v, the integral J Q m (d 3 k/2e 3 )f(0, v) 
can be performed analytically and leads to a function of j3 and ln(fc m /A min ). Its combination with 
the virtual 0(a) R. C. is well defined since, as expected, the I. R. divergences cancel exactly. 

The integration of the second term in Eq. (j5]) is I. R. convergent, so we set A min = (v = 1), 
and replace /(e 7 , v) — /(0, v) —> f(k, 1) — /(0, 1). In this latter expression, the e 7 = k dependence 
is contained in E and, therefore, also in p, (3, and f(k, 1) (cf. Eqs. (0 HJ)). 

Changing the variable of integration from k to E via Eq. (jllh the integral j d 3 k/2k 3 [f(k, 1) — 
/(0, 1)] can be also performed analytically. In particular, we note that the result of this integration 
includes the integral 

1=1 ^ln(l-t)/(t 2 -m 2 lE 2 fl 2 , (6) 

J m/E 
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for which we have derived the compact analytic representation 

/ = -Li 2 (2/3/(1 + $)\ + tanh" 1 /? In (2/3 2 /(l + $)) , (7) 

where z 

U 2 {z) = - [ (dt/t)ln(l-t) (8) 
Jo 

is the dilogarithm. 

As mentioned before, the d 3 k integration of the two remaining contributions, namely the 
interference of the two terms in the I. B. amplitude and the absolute value square of the second, 
are I.R. convergent. Changing again the variable of integration from k to E via Eq. (jlj), they can 
also be performed analytically. 

Combining all the 0(a) contributions discussed above, we obtain the radiatively corrected 
expression for the v e spectrum: 



dP v = dP° 

where 

is the zeroth order spectrum 
h(E,E ) 

tanh _1 /3 



V m / 4 



l + (a/2n)h(E,E ) 
dP° = ApEF(Z,E) K 2 dK 

I tanh _1 /3 



(9) 
(10) 
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2E(3 



m 



2tanh _1 /3 
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is the proton mass, and E,p, f3,Li 2 (z) are defined in Eqs. (!2|8|) . In Eq. (|T0|) . F(Z,E) is the 
Fermi Coulomb function evaluated at E = E and A is a constant independent of K. In neutron 
/9-decay, for example, A = (gy + 3g\)/2n 3 , where gy = GpV u d and g& are the vector and axial 
vector coupling constants. In analogy with the e~ spectrum case PQ, Eqs. (!9"] fTT|) can be applied to 
allowed nuclear (3 decay using an independent particle model of the nucleus. In that case g v + 3g\ 
is replaced by gy\Mp\ 2 + g 2 A \McT\ 2 ■, where Mp and Mqt are the Fermi and Gamow- Teller matrix 
elements. In analogy with g(E,E ), Eq. f llip is valid for both v e and v e transitions. 

Aside from the corrections (a/2n)g(E, E Q ) and (a/2n)h(E, E ), in the Standard Theory there 
are large additional electroweak corrections which give very important contributions to the Ft- 
values of the superallowed Fermi transitions and the neutron lifetime. However, they are constants 
independent of K and E, so they don't affect the shape of the v e and e~ spectra. For this reason, 
they are not included in the present discussion. 

Eq. f lTT]) simplifies considerably in the m — > — > 1) limit. Noting that tanh _1 /3 = 



In 



£■(1 + (3)/m and Li 2 (l) = 7r 2 /6, we readily find: 



lim h(E, E ) 

m— >0 



3 In 



m p 

2Er 



23 

T 



4vr 2 



31n(l-y), 



(12) 
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where y = K/Eq. Thus, the m — > limit of h(E,Eo) is convergent and has a very simple form. 
This is in sharp contrast with the behavior of g(E,E ), the function that describes the R. C. to 
the e~(e + ) spectrum. In fact, in the relativistic m 2 « E 2 limit 



g{E,E ) -> 3 In 



2E 
+4 (In x - 1 
2E 



3 _ 2vr2 

4 ~ IT 

1 — x 3 , f 1— x 
h In 

3x 2 V x 



(1-x) 2 , 
+ ^ In x 



+ ln 



m 



4(1 - x) 
3x 



-3 



(1 - x)< 
6x 2 



+ 4 In 



6x 2 
1 — x 

x 



(13) 



where a; = E/ Eq. Eq. f fl3|) can be gleaned from the expression between curly brackets in Eq. (4.4) 
of Ref. [3], after subtracting 61n(A/m p ) + 9/4 (cf. footnote 23 in Ref. PQ). Thus, the correction 
to the electron spectrum diverges logarithmically in the m — > limit! On the other hand, if 
one multiplies Eq. ( 1T3|) by x 2 (l — x) 2 , which describes the energy dependence of the zeroth order 
probability in the m — > limit, and integrates over < x < 1, the cofactor of ln(2£ , / m ) 
vanishes, so that the correction to the inverse lifetime is indeed free from the electron "mass 
singularity" [3|12|I13] . 

The fact that h(E, E ) is convergent in the m — > limit while g(E, E ) is not, can be explained 
by the following theoretical argument. For given K, in the m — > limit all collinear e — 7 
configurations become energy degenerate with energy Eq — K and generally give rise to mass 
singularities (cf. Eq. fl3])). An elementary, but powerful theorem in Quantum Mechanics leads to 
the conclusion that these singularities are cancelled in the power series expansions of transition 
probabilities if the latter are summed over an appropriate ensemble of such degenerate states [13] . 
In deriving the corrections to the z/ e (z/ e ) spectrum, we perform the d 3 p and d 3 k integrations, so 
indeed we sum over the set of collinear e — 7 configurations that become energy degenerate in 
the m — > limit. Thus, according to this theorem, h(E,E Q ) should be free of lnm singularities, 
as we find. In contrast, this is not the case in the derivation of the electron spectrum, since the 
d 3 p integration is not carried out. As a consequence, g(E,E ) (Eq. ( TT3"|) ) explicitly exhibits a 
lnm singularity. Analogous examples of mass singularities in the differential spectrum, and their 
cancellation in the lifetime, integrated asymmetry and some partial decay rates in muon decay are 
extensively discussed in Ref. [3]. 

A useful consistency check of Eqs. f lT2|T3|) is to multiply Eq. ffl2|) by y 2 (l — y) 2 , Eq. ffT3]) by 
x 2 (l — x) 2 , and integrate over x and y from to 1. The two results should be equal since both lead 
to the correction to the inverse lifetime, a fact that we have verified to high accuracy. Eq. ( TT2l) 
can also be derived directly by neglecting the terms in the integrands that vanish in the m — > 
limit. In that approximation the phase space integrations are very simple, and we find that the 
lnm singularities from the various contributions indeed cancel, leading again to Eq. (fT2|) . 

As an illustrative example, we consider the case E = 8 MeV. Using Eq. (TIT]) , we find that 
the R. C. (a / 27r) h(E, E ) to the v e [y^) spectrum increases the transition probability by 1.17 
% at y = 0.8 (K = 6.4 MeV) and by 0.64 % at y = 0.2 (K = 1.6 MeV). This behavior is 
reflected in the simpler expression in Eq. f JT2|) . which shows a positive, monotonically increasing 
function of y. Using the exact expression for g(E,E ) (Eq. (20b) of Ref. []]), we find that the 
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R. C. (a/2ir)g(E, Eq) to the e (e + ) spectrum decreases the transition probability by 1.13 % at 
x = 0.8 (E = 6.4 MeV) and increases it by 2.90 % at x = 0.2 (E = 1.6 MeV). This behavior can 
be qualitatively understood by noting that the cofactor of In (2Eo/m) in Eq. ( fl3l) is negative for 
x > 0.44 and positive for x < 0.44. As a consequence, the R. C. decreases (increases) the transition 
probability for high (low) energy e~ or e + , an effect that becomes sharper as Eq increases. A similar 
behavior has been found in muon decay [3]. 

In many cases, Eqs. (I12|ll3p provide very good and useful approximations to the exact expres- 
sions for h(E, E ) and g(E, Eq). Their domains of validity are (E — K) 2 » m 2 and E 2 » m 2 , 
respectively. As shown in Ref. pQ, the 0{pi) R. C. to the shape of the e~(e + ) spectrum is not 
affected by strong interactions if terms of 0(aE/M) are neglected. The same conclusion applies 
to the 0(a) R. C. to the v e {ye) spectrum. 



3 Radiative Correction to the Conversion from the e Spectrum to 
the v e Spectrum 



Including the 0(a) R. C, the theoretical expressions for the e and v e spectra in allowed 
/3-decay can be cast in the form 



dPe 

dE 
dP u 

dK 



fe(E,E ) 
U(K,E ) 



where 



fe(E,E ) 

fu(K,E ) 



ApE(E -E) 2 F(Z,E) 1 + —g(E } E ) 

In 



a 



ApE K F(Z, E) l + —h(E,E ) 

2tc 



a 



(14) 
(15) 

(16) 
(17) 



and E and p are defined in Eq. (|2]). Comparing Eqs. ( TT6lfT7l) and recalling that E = Eq — K, we 
see that 



f u (K, Eq) = f e {E, E Q ) 1 + ^h(E, Eq)] / [l + £-g(E, Eq) 



or, neglecting terms of order 0(a 2 



a 



f„(K, Eq) = f e (E, Eq) [1 + — [h(E, Eq) - g(E, Eq] 

Using Eq. (TTTj) and the exact expression for g(E,E ) (Eq. (20b) of Ref. [T]), we find 

h(E, Eq) - g(E, E ) = H _ l L i 2 [JL 

2 p \l + 



(19) 



. tanh 1 B 
+4 I ^— - 1 



P 

tanh _1 /3 



11 + In ( 2 J±\ + 21n/3 + 21n(l - y) - \ny - — 

6 v m J 3 (! -y) 



3-/3 2 



4tanh _1 /3 



6(1 -y)< 



(20) 
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where /3 is defined in Eq. fl2]) and we recall that y = K/E . 

Eqs. fU9|20p describe the conversion from the e~ spectrum in a specific decay into the corre- 
sponding v e spectrum when O(o) R. C. are included. 

As an illustrative example, we consider again the case £"0 = 8 MeV. Using Eqs. ( 1T91I201) we find 

that, for a given e spectrum f e (E, Eq), the R. C. (a/2n) h(E, E ) — g(E, Eq) decreases f u (K, Eq) 

by 1.74 % at y = 0.8 (K = 6.4 MeV) and increases it by 1.77 % at y = 0.2 (K = 1.6 MeV). Thus, 
the correction factor [1 + (a/2ir){h(E, Eq) - g(E, E )) in Eq. (02) is 3.57 % larger at K = 1.6 MeV 
than at K = 6.4 MeV. This behavior reflects the dominant contribution of —g(E,Eo). As a 
consequence, the R. C. decreases (increases) f u (K,E ) for high (low) K, an effect that becomes 
sharper as E increases. 



4 Conclusions 

In summary, we have derived an analytic expression for the 0(a) R. C. to the v e (v e ) spectrum 
in allowed /3-decay (Eqs. (l9lfTTT)). We have found that the m — > limit of this correction is 
convergent and leads to a very simple expression (Eq. ( fi~2l ). This is in sharp contrast to the 
0(a) R. C. to the e~(e + ) spectrum, that diverges as m — > 0, an important difference that we 
have explained on theoretical grounds. We have then used our results to derive the 0(a) R. C. 
to the conversion from the e~ spectrum to the v e spectrum (Eqs. ( IT9|20|) ). a relation that plays 
an important role in reactor studies of neutrino oscillations [9]I10|I11] . A recent study concludes 
that the conversion procedure can be implemented with an error < 1%, provided some conditions 
are met pU]. If this accuracy is achieved, it is likely that the R.C. will play a significant role in 
precision studies of the v e spectrum. 
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